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ABSTRACT 

The  goal  of  this  thesis  is  to  investigate  the  behavior 
of  a  thermoacoustic  prime  mover  filled  with  argon.  Data  was 
collected  for  three  different  operating  conditions:  1)  prior 
to  the  onset  of  acoustic  oscillation,  2)  immediately  after 
onset  and  3)  after  the  prime  mover  had  reached  steady  state, 
at  mean  gas  pressures  of  approximately  80  kPa ,  205  kPa  and  305 
kPa .  Measured  values  of  1/Q  are  compared  to  predictions  of  a 
theory  based  on  Swift's  work  (Swift,  G.  W.,  "Thermoacoustic 
Engines,"  J.  Acoust.  Soc.  Am.  Vol.  84,  October  1988).  After 
onset,  Q  is  determined  from  measurement  of  the  rise  time  of 
the  oscillations.  The  data  set  collected  is  the  first  to  span 
onset  and  results  show  good  agreement  with  theory.  It  is 
concluded  that  the  theory  can  be  used  to  provide  reasonable 
input  for  a  model  to  predict  the  steady  state  waveforms. 
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I.   INTRODUCTION 

A  thermoacoustic  prime  mover  is  a  heat  engine;  it 
converts  stored  energy  into  useful  work,  in  the  form  of  sound. 
The  basic  operation  and  construction  of  a  prime  mover  can  be 
explained  with  reference  to  Fig.  1.  For  our  purposes,  a  prime 
mover  can  be  considered  to  be  a  closed  resonator  comprised  of 
five  sections.  A  temperature  gradient  is  applied  across  the 
prime  mover  stack  via  the  ambient  and  hot  heat  exchangers. 
(The  terms  "ambient"  and  "hot"  refer  to  the  temperatures  of 
the  heat  exchangers  and  the  portions  of  the  resonator  adjacent 
to  them.)  The  (prime  mover)  stack  and  heat  exchangers  are 
parallel  plates  spaced  by  a  few  thermal  penetration  depths  and 
oriented  such  that  the  plates  lie  along  the  axis  of  the 
resonator.  The  hot  and  ambient  sections  are  circular  tubes  -- 
open  at  the  heat  exchanger  end  and  closed  at  the  other  end. 
The  net  absorption  coefficient  of  the  prime  mover  is  the  sum 
of  the  thermal  and  viscous  wall  losses  in  each  of  the  five 
sections.   Attenuation  in  the  gas  is  negligible. 

As  the  temperature  difference  imposed  across  the  stack 
increases  from  zero,  the  thermal  losses  in  the  stack  decrease 
and  eventually  become  negative,  representing  gain.  At  a 
sufficiently  high  temperature  difference,  the  gain  of  the 
stack  overcomes  the  remaining  losses  and  the  prime  mover 
spontaneously  goes  into  self  oscillation.   This  phenomenon  is 
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known  as  onset.  Once  onset  is  reached,  the  acoustic  amplitude 
of  the  resonator  grows  rapidly,  typically  reaching  a  few 
percent  of  the  mean  pressure.  This  discussion  leads  to  the 
conclusion  that  three  distinct  regions  of  prime  mover  behavior 
are  of  interest:  1)  below  onset,  2)  above  onset  during  the 
buildup  of  oscillations  to  steady  state,  and  3)  steady  state. 
The  purpose  of  this  thesis  is  to  investigate  each  of  these 
regions  in  an  argon-filled  prime  mover. 

This  thesis  is  one  in  a  sequence  that  have  investigated 
thermoacoustic  prime  movers.  The  goal  of  this  research  is  to 
understand  the  behavior  of  prime  movers  in  each  of  the  three 
regions  of  operation.  In  1989,  Lin  began  a  quantitative  study 
of  a  helium-filled  prime  mover.  He  investigated  the  region 
below  onset  as  well  as  documenting  the  nature  of  steady  state 
oscillations.  Lin  found  that  the  steady  state  oscillations  in 
our  prime  mover  are  highly  nonlinear,  exhibiting  a  rich 
harmonic  structure.  (Lin,  1989)  Therefore,  it  became  clear 
that,  to  completely  understand  the  steady  state  region,  it  is 
also  necessary  to  understand  the  generation  of  harmonics  and 
energy  dissipation  in  high  amplitude  standing  waves.  In  1991, 
Chen  began  such  an  investigation  in  a  simple  standing  wave 
tube  (Chen,  1991).  Also  in  1991,  Bowers  investigated  the 
initial  buildup  of  oscillations  above  onset  (Bowers,  1991). 

As  stated  earlier,  the  purpose  of  this  thesis  is  to 
investigate  the  three  regions  of  behavior  in  an  argon-filled 
prime   mover.     One   reason   for   this   is   to   add   to   our 


understanding  of  prime  movers  by  using  a  gas  other  than 
helium.  Also,  Lin  and  Bowers  used  slightly  different  prime 
movers  so  direct  comparison  of  their  data  is  not  possible.  We 
will  use  the  same  prime  mover  configuration  below  and  above 
onset  so  that  we  can  have  a  complete  data  set  from  zero 
temperature  difference  through  onset  and  beyond. 

Following  a  discussion  of  applicable  theory,  the 
experimental  setup  and  procedures  used  to  collect  data  will  be 
presented.  The  results  and  their  analysis  follow,  along  with 
a  discussion  of  sources  of  error.  Finally,  a  summary  of  this 
experiment  and  its  conclusions  will  be  made  along  with 
recommendations  for  subsequent  research  in  this  area. 


II.   THEORY 

Theoretically,  the  performance  of  a  prime  mover  can  be 
investigated  through  the  calculation  of  its  acoustic  power 
output.  However,  below  onset  and  during  the  initial  buildup 
of  oscillations,  it  is  easier  to  measure  the  quality  factor  of 
the  prime  mover. 

The  net  acoustic  power  output,  W,  can  be  related  to  the 
quality  factor,  Q.  Below  onset,  Q  is  easily  determined  from 
the  prime  mover's  frequency  response.  Just  above  onset,  Q  can 
be  determined  from  the  rise  time  of  the  oscillations.  This 
section  contains  a  discussion  of  the  theory  necessary  to 
calculate  1/Q  in  both  regions.  This  discussion  is  a  summary 
of  that  found  in  Bowers  (Bowers,  1991,  pp.  10-15)  which  draws 
heavily  from  Swift  (Swift,  1988)  . 

The  reciprocal  of  the  quality  factor,  1/Q,  of  the  prime 
mover  can  be  expressed  as 
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where  a;  is  the  angular  frequency  and  EST,  the  stored  energy  of 
the  resonant  tube.  W  is  determined  by  integrating  the  power 
per  unit  volume,  w(x),  over  the  volume  of  the  prime  mover,  as 
follows, 
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S(x)  is  the  cross-sectional  area  of  the  prime  mover  at  x. 
In  the  stack,  w(x)  is 
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PA  and  k  are  the  peak  pressure  amplitude  and  wave  number  of 
the  standing  wave,  whereas  a  and  a  represent  the  thermodynamic 
sound  speed  and  Prandtl  number  of  the  working  fluid.  j2  =  -1. 
Also,  VTm  is  the  mean  temperature  gradient  along  the  stack 
plates  and  T  is  the  mean  temperature  of  the  plate  and  is  a 
function  of  x.   f  is  defined  as 
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for  parallel  plates.   S  ,     the  viscous  penetration  depth,  is 
defined  as 
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/i,  is  the  fluid's  dynamic  viscosity.   A  similar  relationship 


exists  between   f  ,  6     and  y  ,  although  8     is  defined  as 
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K  is  the  thermal  conductivity  of  the  working  fluid. 

es  defines  the  relative  heat  capacity  of  the  gas  compared 
to  the  plates 
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where  p  and  c  are  the  mean  density  and  isobaric  specific 
heat  of  the  working  fluid  and  ps,  cs  and  Ss  represent  similar 
qualities  for  the  stack  plates.  J  is  one-half  the  stack  plate 
thickness.  In  the  heat  exchangers,  w(x)  is  identical  to 
Equation  (3),  except  ATm  =0.  In  the  ambient  and  hot  ends, 
w(x)  is  given  by  Equation  (3)  with  ATm  =  0  and  1=0  (there 
are  no  plates)  .  Also,  the  functions  f  and  f  take  on 
different  forms  involving  the  Bessel  functions. 

EST  is  determined  by  integrating  the  time  averaged  energy 
density  over  each  prime  mover  section, 
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where  S  is  the  cross-sectional  area  of  the  resonator  and  <p   is 
the  phase  of  the  standing  wave  as  it  enters  each  section. 

These  equations  use  the  properties  of  the  working  fluid 
and  the  initial  conditions  of  the  prime  mover  to  predict  1/Q 
as  a  function  of  VT   across  the  stack. 

m 

The  resonance  frequency  is  calculated  by  neglecting  the 
dependence  of  the  sound  speed,  a,  on  dissipation.  In  other 
words  a  is  assumed  to  be  the  thermodynamic  sound  speed.  Next, 
we  require  that  acoustic  pressure  antinodes  exist  at  each  end 
of  the  prime  mover  which  amounts  to  requiring  that 


\k(x)  dx  -  mi 
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k   depends    on    x   through   k(T  (x)  ) 


III.   EXPERIMENT 

A.   EXPERIMENTAL  SETUP 

The  prime  mover  and  temperature  and  pressure 
monitoring/control  systems  are  same  ones  used  in  Bowers7 
thesis.   Only  pertinent  details  will  be  given  here. 

1.   The  Prime  Mover 

As  discussed  earlier  and  shown  in  Fig.  1,  the  prime 
mover  is  composed  of  five  sections. 

a.  The   Hot   End 

The  hot  end  of  the  prime  mover  consists  of  a  nickel 
tube,  0.63  cm  thick,  5 . 5  cm  long  and  having  an  inner  diameter 
of  3.82  cm.  A  nickel  cap  is  welded  into  the  top  of  this  tube 
to  provide  a  rigid  boundary.  To  monitor  the  hot  end 
temperature,  a  small  hole  is  drilled  into  the  end  cap  and 
fitted  with  an  access  tube  for  a  K-type  thermocouple.  The 
lower  end  of  this  assembly  holds  the  hot  heat  exchanger. 

b.  The   Hot   Heat   Exchanger 

The  hot  heat  exchanger  is  also  composed  of  nickel 
and  like  the  hot  end,  has  an  inner  diameter  of  3.82  cm.  25 
nickel  plates,  0.045  cm  thick  and  0.762  cm  long  are  mounted 
along  the  heat  exchanger's  longitudinal  axis.  The  spacing 
between  these  plates  is  0.104  cm.  304  stainless  steel  spacers 
mounted  to  one  side  of  each  plate  help  maintain  this  spacing 
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and  prevent  deformation  under  high  heat  conditions.  The 
hot  end/heat  exchanger  assembly  is  welded  to  the  prime  mover 
stack . 

c.  The   Prime   Mover  Stack 

The  prime  mover  stack  is  housed  in  a  thin  walled 
304  stainless  steel  tube,  3 . 5  cm  long,  with  an  inner  diameter 
of  3.82  cm  and  a  wall  thickness  of  0.05  cm.  The  stack 
consists  of  35,  304  stainless  steel  plates,  0.028  cm  thick, 
spaced  0.077  cm  apart.  As  in  the  hot  heat  exchanger,  the 
spacing  between  these  plates  is  maintained  by  using  304 
stainless  steel  spacers  mounted  to  one  side  of  each  plate. 

d .  The   Ambient   Heat   Exchanger 

A  constant  ambient  temperature  is  maintained  at  the 
bottom  of  the  stack  by  using  2  identical  copper  heat 
exchangers.  Each  is  1.02  cm  long  and  consists  of  25  copper 
plates,  0.025  cm  thick  and  spaced  0.105  cm  apart.  The  inner 
diameter  of  each  heat  exchanger  is  3.82  cm.  The  two  heat 
exchangers  are  spaced  0.15  cm  apart.  The  top  heat  exchanger 
is  mounted  in  contact  with  the  lower  end  of  the  prime  mover 
and  positioned  so  that  its  plates  form  right  angles  with  the 
stack  plates.  This  maximizes  fluid  flow  out  of  the  stack  by 
ensuring  that  the  heat  exchanger  plates  do  not  block  the  gaps 
between  the  prime  mover  plates.  The  second  heat  exchanger  is 
not  necessary  for  the  purposes  of  this  thesis.  It  will  be 
used  in  later  investigations. 
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A  constant  working  temperature  is  maintained  within 
the  heat  exchangers  by  mounting  them  in  a  brass  water  jacket. 
This  water  jacket  also  contains  a  small  hole  used  for 
pressurizing  the  prime  mover.  The  hot  end  heat 
exchanger/stack  assembly  is  soldered  to  the  top  of  the  water 
jacket. 

e.  The  Ambient   End 

The  ambient  end  of  the  prime  mover  consists  of  a 
3.82  cm  ID,  88.11  cm  long  copper  tube,  flanged  at  the  bottom 
to  accommodate  an  end  cap.  Two  different  end  caps  are  used. 
Below  onset  the  end  cap  houses  an  electret  driver  and  an 
electret  microphone,  so  that  the  frequency  response  of  the 
tube  can  be  determined.  The  end  cap  used  to  collect  data 
after  onset  holds  only  a  microphone  to  record  the  amplitude  of 
the  standing  wave. 

A  butterfly  valve  is  located  38.0  cm  from  the 
bottom  of  the  ambient  heat  exchanger.  Closing  the  butterfly 
valve  effectively  shortens  the  length  of  the  prime  mover  and 
requires  a  greater  differential  temperature  to  drive  the  prime 
mover  to  onset.  With  the  butterfly  valve  closed,  a 
differential  temperature  across  the  stack  sufficient  for  onset 
can  be  established  without  actually  generating  a  standing  wave 
(Bowers,  1991,  p.  18). 
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2.   Control  and  Monitoring  Systems 
a.  Temperature 

The  hot  end  of  the  prime  mover  is  fitted  with  an 
Omega  Engineering  HBA  Model  202040  heater.  Power  is  supplied 
to  the  heater  by  a  General  Radio  Type  100-Q  variable 
transformer  which  controls  the  heater  voltage.  The  heater 
coils  are  wrapped  tightly  around  the  hot  end  of  the  prime 
mover  and  stainless  steel  shim  stock  is  used  to  maintain  good 
thermal  contact.  A  Neslab  RTE-110  constant  temperature  bath 
supplies  cooling  water  to  the  water  jacket  and  to  plastic 
tubing  which  encircles  the  entire  length  of  the  ambient  end  of 
the  prime  mover. 

Temperatures  along  the  length  of  the  prime  mover  are 
monitored  by  a  Keithley  Model  740  System  Scanning  Thermometer. 
A  K-type  thermocouple  probe  is  inserted  into  the  access  hole 
in  the  hot  end  until  that  it  makes  contact  with  the  hot  heat 
exchanger.  Three  E-type  thermocouple  junctions  are  affixed  to 
the  outside  of  the  ambient  end  of  the  prime  mover  with 
cyanoacrylate  adhesive.  These  junctions  are  placed  below  the 
ambient  heat  exchanger,  below  the  butterfly  valve  and  above 
the  pressure  transducer  flange  at  the  bottom  of  the  ambient 
end. 
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jb .  Pressure 

The  prime  mover  in  this  thesis  uses  argon  gas  as 
the  working  fluid.  The  gas  pressure  is  controlled  by  a 
Matheson  8H-580  pressure  regulator  with  a  working  range  of 
0  -  1400  kPa .  Prior  to  pressurization,  the  prime  mover  is 
evacuated  using  a  Welch  Model  1400  Duo-Seal  Vacuum  Pump  for 
several  hours  to  ensure  all  residual  gases  are  removed. 
Pressure  within  the  prime  mover  is  monitored  via  the  pressure 
sensing  line  in  the  waterjacket  assembly.  An  external  tubing 
and  valve  network  allows  the  pressure  sensing  line  to  be 
isolated  from  the  gas  supply  cylinder  and  vacuum  pump.  The 
pressure  is  monitored  by  an  Omega  PX304-300AV  pressure 
transducer  having  a  sensitivity  of  3PSIA/MVDC.  The  output  of 
the  pressure  transducer  is  read  by  a  Hewlett  Packard  3457A 
digital  multimeter. 

B.   EXPERIMENTAL  PROCEDURE 

As  power  is  supplied  to  the  hot  end  of  the  prime  mover, 
the  differential  temperature,  ATm,  increases  from  zero  to  a 
level  sufficient  to  bring  about  onset.  Data  was  acquired  in 
the  three  regions  of  prime  mover  operation  discussed  in  the 
Introduction. 

1.   Below  Onset 

a.  Data   Acquisition 

A  schematic  diagram  of  the  apparatus  used  to 
measure  frequency  response  of  the  prime  mover  below  onset  is 
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shown  in  Fig.  2.  Data  acquisition  is  controlled  using  GPIB 
hardware  and  a  program  written  in  Microsoft  QuickBasic  and 
executed  on  a  Standard  286  PC-AT.  After  entering  the  expected 
resonance  frequency,  desired  bandwidth,  number  of  frequency 
increments  and  starting  system  pressure,  the  PC  directs  a 
Hewlett-Packard  3324A  Synthesized  Frequency  Generator  to 
output  a  10  Vpp  sine  wave  beginning  at  the  lower  end  of  the 
desired  bandwidth.  This  signal  is  sent  to  the  reference  input 
of  a  Stanford  Research  Systems  530  Lock-In  Amplifier,  a 
Kikusui  COS6100A  Oscilloscope  and  a  Techron  7520  Power  Supply 
Amplifier  which  amplifies  the  signal  and  sends  it  to  the 
driver . 

The  driver  is  a  locally  constructed  electret  type 
which  sets  up  a  standing  wave  in  the  prime  mover.  An  electret 
microphone  is  used  to  sense  the  amplitude  of  the  standing 
wave.  The  output  of  the  microphone  is  first  directed  through 
a  locally  constructed  preamplifier  which  boosts  the  signal  by 
a  factor  of  100.  The  amplified  signal  is  then  sent  to  the 
lock-in  amplifier  and  oscilloscope. 

The  oscilloscope  provides  a  visual  display  of  the 
reference  signal  and  the  received  signal,  while  the  lock-in 
amplifier  filters  the  received  signal  in  a  narrow  bandwidth 
around  the  reference  frequency.  It  sends  the  in-phase  and 
quadrature  components  of  the  signal  to  the  HP  3457A 
multimeter.  These  values  are  read  by  the  PC,  via  a  GPIB 
interface,  and  the  amplitude  of  the  signal  is  determined. 
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After  the  acoustic  amplitude  data  and  the 
temperature  of  each  thermocouple  are  written  to  a  worksheet, 
as  shown  in  Fig.  3,  the  frequency  is  incremented.  This 
procedure  is  repeated  until  the  end  of  the  bandwidth  is 
reached,  at  which  time  the  ending  system  pressure  is  recorded. 
Data  was  recorded  for  both  the  fundamental  and  second  modes  of 
the  prime  mover  at  the  three  mean  gas  pressures  of  interest. 
b.    Data   Processing 

Q  is  obtained  by  applying  the  same  fitting  function 
Lin  used  to  determine  the  resonant  response  of  a  simple  tube. 
That  function  is 


P  = 


J_ 


cot {kL) 


kL 


(9) 


where  P  is  the  normalized  amplitude  of  the  acoustic  pressure 
in  the  tube  and  A  is  a  scaling  constant.  The  subscript  "o" 
indicates  a  measured  resonant  value.   kL  is  defined  as 


kL  -  n 


w 


J 


(20) 


(10) 


By  performing  a  least  squares  fit  of  P  versus  f,  values  for 
f  ,  A  and  Q  can  be  determined.  (Lin,  1989,  p.  22)  Data  is 
processed  using  macro  routines  in  Microsoft's  Excel  program. 
Figure  4  shows  a  typical  plot  resulting  from  the  fitting 
procedure . 
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Figure  4.   Sample  plot  of  data  obtained  prior  to  onset. 

To  compensate  for  any  loss  of  pressure  during  the 
run,  the  mean  gas  pressure  is  obtained  by  averaging  the 
starting  and  ending  system  pressures.  The  mean  hot  and 
ambient  end  temperatures  are  the  averages  of  the  temperatures 
recorded  at  each  incremental  freguency. 
2 .   Determining  Onset 

Unlike  helium,  which  exhibits  a  well-defined  onset 
temperature  when  used  as  the  working  fluid,   it  was  more 
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difficult  to  identify  the  experimental  onset  temperature  in 
argon . 

As  onset  was  approached,  a  very  low  amplitude  acoustic 
wave  was  detected  by  the  sensing  equipment.  The  amplitude  of 
this  wave  increased  slowly,  requiring,  in  some  cases,  several 
minutes  to  build  to  steady  state.  The  exact  reason  for  this 
behavior  is  unknown.  The  temperature  differences  required  for 
onset  are  shown  in  Table  1.  These  values  were  determined  by 
slowly  increasing  ATm  until  onset  occurred,  then  observing  the 
behavior  of  the  prime  mover  at  several  temperatures  around 
this  reference  point.  AT0nset  is  the  mean  value  of  the 
temperatures  for  which  onset  was  observed.  At  onset,  1/Q  =  0. 


Pressure    (kPa) 

80 

205 

305 

AT0nset         CC) 

366.4 

381.7 

403  .5 

Table  1.   SUMMARY  OF  ONSET  TEMPERATURES 
3.   Above  Onset 

a.  Data   Acquisition 

A  schematic  diagram  of  the  apparatus  used  to 
collect  data  above  onset  is  shown  in  Fig.  5.  The  pressure 
transducer  used  for  this  phase  of  data  collection  is  an 
ENDEVCO  Model  8510B-5  piezoresistive  microphone  with  a  linear 
response  range  of  +  5  PSI.  Data  is  collected  above  onset 
using  the  procedures  outlined  in  Bowers'  thesis.   With  the 
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is  established  and  allowed  to  stabilize.  When  the  butterfly 
valve  is  opened,  acoustic  oscillations  begin  to  build  up  in 
the  prime  mover.  The  resulting  acoustic  wave  is  sensed  by  the 
microphone  and  routed  through  a  locally  constructed, 
differential  amplifier  having  a  gain  of  10.  This  signal  is 
sent  to  a  Nicolet  310  digital  oscilloscope  for  recording  and 
storage.  (Bowers,  1991,  p  20)  The  oscilloscope  is  triggered 
to  begin  recording  data  immediately  after  onset  occurs, 
ensuring  that  the  acoustic  signal  is  small.  This  minimizes 
the  presence  of  harmonics  in  the  recorded  signal  and  ensures 
that  the  temperature  distribution  in  the  stack  is  unaltered  by 
significant  acoustic  heat  transport. 
b.    Data   Processing 

The  Q  of  the  prime  mover  above  onset  is  determined 
by  the  relationship, 

0-^,  (ID 

where  go  is  the  angular  resonance  frequency  and  fi  is  the 
temporal  absorption  coefficient.  To  ascertain  the  value  of 
these  variables,  the  data  collected  above  onset  is  digitally 
processed  using  the  techniques  developed  by  Bowers  (Bowers, 
1991,  pp  22-27)  . 

The  binary  signal  stored  by  the  Nicolet  310  is 
converted  to  ASCII  code  to  be  processed  using  Mathworks  MATLAB 
program.    Initially,  a  discrete  Fourier  transform,  DFT,  is 
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performed  and  the  positions  of  the  resonance  frequency  and  any 
harmonics  are  determined.  A  raw  waveform  is  shown  in  Fig.  6, 
while  Fig.  7  shows  the  DFT  of  this  data. 

The  original  signal  is  multiplied  by  a  complex 
exponential  term  which  effectively  shifts  the  fundamental 
frequency  to  DC  (i.e.  0  Hz).  A  5th  order,  50  Hz  Butterworth 
filter  is  then  applied  to  the  shifted  data  set  to  isolate  the 
fundamental  frequency.  Figure  8  shows  a  plot  of  the  shifted, 
filtered  data  revealing  that  only  the  envelope  of  the  original 
data  remains.  As  shown  in  Bowers'  thesis,  the  end  effects  of 
the  Butterworth  filter  can  be  avoided  by  disregarding  the 
first  300  and  last  300  points  of  filtered  data.  The 
exponential  envelope  of  the  original  waveform  is  linearized  by 
taking  the  logarithm  of  the  normalized  amplitude,  so  that 
aleast  mean  squares  fit  can  be  performed.  The  slope 
determined  by  this  fit  corresponds  to  the  desired  quantity,  (3. 
4.   Steady  State 

a.  Data   Acquisition 

Collection  of  data  once  the  prime  mover  has  reached 
steady  state  uses  the  same  experimental  setup  as  was  used  to 
collect  data  above  onset.   With  the  butterfly  valve  open  and 
a  sufficient  AT  established  to  induce  self  oscillation,  the 

prime  mover  is  allowed  to  reach  equilibrium  and  the  resulting 
waveform  is  recorded  by  the  Nicolet  oscilloscope. 
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Figure  6.   Sample  plot  of  data  collected  during  onset 
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In  order  to  prevent  aliasing,  the  digital  sampling 
rate  must  conform  to  the  Nyquist  criteria, 


£ Ny guise    ^  ^ -^max  ' 


(12) 


where  fN  ist  is  the  minimum  sampling  frequency  and  fmax  is  the 
highest  frequency  of  interest  in  the  signal  being  sampled. 
(Strumm,  1988,  p  54)  Since  the  fundamental  frequency  of  the 
prime  mover  is  approximately  160  Hz  and  there  will  be  at  most 
10  harmonics  detectable  in  the  steady  state,  f  =  1760  Hz  and 
f„   .  .  =  3520  Hz.   This  corresponds  to  a  minimum  sampling  rate 

Nyquist  r  r 

of    308    /lis    per   point.       We    used    200    /is    per    point. 
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Since  we  are  interested  in  the  absolute  pressure  amplitude 
of  the  harmonics  of  the  steady  state  waveform,  a  known  test 
signal  is  also  recorded,  using  the  same  oscilloscope  settings. 
The  DFT  of  this  signal  will  provide  a  reference  when  analyzing 
the  original  data. 

b.    Data  Processing 

Using  Mathworks  MATLAB  program,  a  DFT  of  the 
reference  signal  and  the  recorded  waveform  are  performed,  as 
shown  in  Figs  9  and  10.  In  the  interest  of  clarity,  Fig.  9 
shows  only  the  first  1000  points  of  the  original  signal's  4000 
point  DFT,  while  only  the  first  half  of  the  DFT  of  the 
recorded  waveform  is  shown  in  Fig.  10. 
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Figure  9.   Plot   of   DFT  of   reference   signal   used   for 
analysis  of  steady  state  data. 
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IV.   RESULTS 

Data  were  taken  at  80  kPa ,  205  kPa  and  305  kPa  for  AT  's 
ranging  from  O'C  to  approximately  425°C,  the  self-imposed  safe 
working  limit  for  the  prime  mover  assembly.  Below  onset,  data 
for  the  fundamental  and  second  modes  were  collected  at  all 
three  pressures.  However,  due  to  temperature  constraints, 
only  the  two  lower  pressures  were  investigated  above  onset  and 
at  steady  state. 

1.   Results  of  1/Q  Measurements 

Figures  11  and  12  show  the  combined  results  of  the 
data  collected  below  onset  and  immediately  above  onset  for 
pressures  of  80  kPa  and  205  kPa .  Fig.  13  shows  the  results 
obtained  below  onset  at  305  kPa .  For  discussion  purposes,  it 
should  be  pointed  out  that  1/Q  is  proportional  to  the  net 
attenuation  in  the  prime  mover. 

The  first  point  of  discussion  is  that  the  data  in 
Figs.  11  and  12  show  a  smooth  transition  from  below  onset,  up 
to  and  beyond  onset.  This  feature  indicates  that  the  methods 
of  determining  1/Q  in  each  region  are  consistent.  These  data 
are  the  first  to  show  the  progression  from  below  onset  to 
above  onset.  It  should  be  pointed  out  that  the  above  onset 
data  were  obtained  at  low  acoustic  amplitudes.   This  ensures 
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Figure  11.  Plot  of  combined  1/Q  vs  ATm  data  obtained  below 
and  immediately  after  onset.  P  =  80kPa. 
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Figure  12.   Plot  of  combined  1/Q  vs  ATm  data  obtained  below 
oand  immediately  after  onset.  P  =  205  kPa . 
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Figure  13.  Plot  of  1/Q  vs  ATm  for  data  obtained  below  onset 
in  the  fundamental  mode.  P  =  305  kPa . 


that  the  effects  of  acoustic  heat  transport  can  be  ignored, 
just  as  they  were  below  onset.  The  results  for  the  second 
mode  are  shown  in  Figs.  14,  15  and  16.  All  of  these  data  were 
obtained  below  onset. 

The  agreement  between  data  and  theory  is  good, 
however,  there  is  a  tendency  to  underpredict  the  attenuation. 
This  trend  is  especially  evident  in  the  second  mode.  One 
possible  explanation  for  this  is  the  fact  that  the  model  does 
not  include  the  presence  of  the  butterfly  valve.  The 
magnitude  of  this  effect  should  be  roughly  proportional  to  the 
additional  surface  area  introduced  by  the  valve.  The  valve  is 
a  circular  disk  of  negligible  thickness  with  a  diameter  of 
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Figure  14.   Plot  of  1/Q  vs  ATm  for  data  obtained  below  onset 
in  the  second  mode.   P  =   80  kPa . 
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Figure  15.  Plot  of  1/Q  vs  ATm  for  data  obtained  below  onset 
in  the  second  mode.  P  =  205  kPa . 
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Figure  16.  Plot  of  1/Q  vs  VTm  for  data  obtained  below  onset 
in  the  second  mode.  P  =  305  kPa . 


approximately  38  cm.  The  total  surface  area  is,  therefore, 
approximately  24  cm2.  The  total  surface  area  of  the 
primemover,  including  heat  exchangers  and  prime  mover  stack, 
is  approximately  2300  cm2.  Therefore,  the  valve  adds 
approximately  1  percent  more  surface  area  to  the  prime  mover. 
Although  a  more  detailed  calculation  of  the  effect  of  the 
valve  is  needed,  it  is  unlikely  that  it  alone  can  account  for 
all  the  error.  It  should  be  pointed  out  that  the  same  trend 
is  evident  in  helium  data.  (Atchley,  1991) 

Another  possible  explanation  concerns  the  fact  that 
Thot  and  Tcold  are  recorded  at  the  top  of  the  hot  heat  exchanger 
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and  the  bottom  of  the  ambient  heat  exchanger,  respectively. 
These  temperatures  may  not  accurately  reflect  the  actual 
temperature  difference  across  the  stack.   This  suggests  that 

the  actual  AT  is  less  than  that  measured,  an  error  that  would 

in  ' 

in  effect  shift  the  data  to  the  left,  providing  a  much  better 
overall  agreement.  The  magnitude  of  this  error  is  difficult 
to  measure.  Measurements  performed  to  try  and  estimate  the 
error  revealed  that  the  results  were  very  sensitive  to  the 
degree  of  contact  between  the  heat  exchangers  and  the  stack. 
A  discrepancy  of  +  10°C  is  certainly  not  unreasonable. 
However,  without  knowing  the  degree  of  contact  in  the  prime 
mover,  any  estimates  of  the  magnitude  of  the  error  are 
speculation. 

Another  point  of  comparison  has  to  do  with  the 
dependence  of  the  difference  between  theory  and  experiment  on 
ATm.  Figures  11  and  12,  and,  to  a  lesser  extent,  Fig.  13, 
show  that  the  agreement  is  worse  at  higher  values  of  AT  . 
This  trend  is  absent  in  Figs.  14,  15  and  16.  A  possible 
explanation  for  this  trend  may  the  irregularities  in  stack 
geometry.  The  theory  assumes  plate  spacing  is  uniform, 
however,  this  is  not  the  case.  One  would  expect  this  effect 
to  be  most  pronounced  when  6  is  comparable  to  y  .  Recalling 
the  definition  of  S  from  Eqn .  (6),  this  would  be  true  for  low 
pressures  and  frequencies.  This  follows  exactly  the  trend 
noted. 
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2 .   Other  Sources  of  Error 

In  addition  to  the  specific  sources  of  error 
previously  mentioned,  there  are  others  associated  with  the 
apparatus  and  procedures  which  apply  to  all  data  collected. 

a.    Monitoring  Systems 

The  measurements  of  mean  gas  pressure,  P  ,  and 
system  temperatures  are  subject  to  instrumentation  errors  and 
errors  due  to  averaging.  T  and  Th  are  used  to  determine  ATm, 
while  T  and  P  are  inputs  to  the  equation  used  for  theoretical 

c        m  c  ^ 

prediction  of  prime  mover  response. 

For  example,  when  investigating  the  response  of  the 
fundamental  mode  at  P  w  205  kPa  and  AT  ~  150°C,  the  average 
value  of  T.  over  the  course  of  data  collection  is  172.2  + 

n  — 

0.5  °C  and  T   =  20.2  +  0.03  °C   As  a  result  AT  =  152.0  + 

c  —  m  — 

0.5  °C,  an  error  of  0.33  percent.  P  is  the  average  of  the 
mean  system  pressures  at  the  start  and  finish  of  each  data 
collection  run.  Actual  system  pressure  varies  over  the  course 
of  collecting  data  due  to  small  leaks  in  the  prime  mover. 
Typically,  the  system  lost  one  percent  of  pressure  at  205  kPa, 
two  percent  at  305  kPa  and  gained  0.2  percent  at  75  kPa . 

The  prediction  of  prime  mover  behavior  for  a  given 
pressure  and  mode  requires  T  and  P  as  input.  These 
parameters  are  obtained  by  averaging  their  values  for  an 
entire  data  set.   For  example,  in  the  fundamental  mode  at 
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approximately  205  kPa ,  Tc  =  20.3  +  0.3°C  and  Pm  =  203.8  +  2.1 
kPa,  differences  of  1.5  and  1.0  percent,  respectively. 

Furthermore,  the  sensors  used  to  monitor  these 
parameters  are  subject  to  design  errors.  According  to  the 
manufacturer's  specifications,  the  error  in  K-type 
thermocouples  is  +  5°C  while  the  E-type  thermocouples  are 
subject  to  an  error  of  +  2°C.  The  Omega  pressure  transducer 
has  an  error  of  +  0.5  percent,  based  on  manufacturer's  data. 

b.    Determining   the   Resonance   Frequency 

Prior  to  onset,  the  fitting  function  used  to 
determine  the  resonance  frequency,  f  is  accurate  to 
approximately  0.1  percent.  Above  onset,  f  is  determined  by 
analyzing  the  DFT  of  the  sampled  waveform.  A  sampling  rate  of 
5  jiisec  for  a  4000  point  record  corresponds  to  a  frequency 
range  of  2000  Hz.  This  results  in  an  accuracy  of  0.5  Hz  per 
point  or  0.3  percent  for  f  =  160  Hz.  A  larger  error  in  the 
frequency  occurs  as  a  result  of  calculating  f  .  It  was  stated 
in  Chapter  II  that  the  computation  neglects  dispersion  due  to 
attenuation.  The  result  is  that  the  theory  slightly 
overpredicts  fQ.  However,  the  maximum  error  is  approximately 
2  percent. 

3.   Steady  State 

Figures  17  through  20  show  the  waveform  and  DFT  of 
data  collected  after  the  prime  mover  had  reached  steady  state 
above  onset  at  approximately  80  kPa .    Figs.   17   and   18 
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correspond  to  a  ATm  of  382.8  °C,  while  Figs  19  and  20 
correspond  to  a  ATm  of  419. 8°C.  An  extensive  set  of  data  is 
presented  in  the  Appendix.  It  is  seen  from  these  figures  that 
the  waveforms  become  highly  nonlinear  at  higher  temperature 
differences.  These  results  are  consistent  with  Lin's  data 
obtained  in  helium.  A  future  goal  of  this  research  is  to 
predict  these  waveforms. 
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Figure  18.  Plot  of  DFT  of  data  obtained  at  steady  state 
P  =  8  0.1  kPa.  AT  =  3  8  2.8  °C. 
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V.   CONCLUSIONS 

The  purpose  of  this  thesis  was  to  investigate  the  behavior 
of  an  argon-filled  thermoacoustic  prime  mover.  The  reasons 
for  this  investigation  were  two-fold.  First,  the  data 
collected  below  and  immediately  above  thermoacoustic  onset 
allowed  a  verification  of  theoretical  predictions  begun  in 
previous  research.  Since  the  previous  work  used  helium  as  the 
working  fluid,  this  experiment  varied  the  input  parameters  by 
using  argon,  instead.  Additionally,  this  data  is  the  first  to 
span  onset.  Our  results  indicate  that  the  theory  provides  a 
good  agreement  with  measurement,  although  some  discrepancies 
do  exist.  At  lower  mean  gas  pressures,  the  prime  mover's 
behavior  begins  to  deviate  from  theory  at  higher  values  of 
ATm.  This  deviation  may  be  due  to  imperfections  in  the  prime 
mover  stack  geometry.  There  is  also  a  tendency  to 
underpredict  the  attenuation. 

Data  was  also  obtained  at  steady  state.  These  data 
show  that  the  waveforms  are  highly  nonlinear.  One  of  the 
ultimate  goals  of  our  research  is  to  be  able  to  predict  the 
shape  of  these  waveforms.  At  this  time  we  are  working  with 
a  model  that  predicts  the  waveforms  of  high  amplitude  standing 
waves  in  simple  tubes.  The  inputs  to  the  model  are  the 
resonance  frequency  and  the  Q  of  the  modes  of  the  tube.  Once 
a  prime  mover  has  gone  into  oscillation,   it  is  no  longer 
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possible  to  measure  these  parameters.  We  will  have  to  rely  on 
theoretical  values. 

The  good  agreement  between  the  measured  and  predicted 
values  of  1/Q  indicate  that  we  will  be  able  to  provide  a 
reasonably  good  input  for  the  model.  This  discussion  leads  to 
a  somewhat  subtle,  but  important,  point.  Below  onset,  Q  is 
determined  by  fitting  the  measured  frequency  response  to  the 
ideal  frequency  response  of  a  simple  tube.  The  quality  of  the 
agreement  indicates  that,  to  a  good  approximation,  the  prime 
mover  behaves  like  a  simple  tube  having  a  uniform  cross 
section  and  temperature  distribution! 
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APPENDIX 
PLOTS  OF  STEADY  STATE  DATA 
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Figure  22.  Plot  of  DFT  of  data  obtained  at  steady  state 
P  =  81.1  kPa.  AT  =  398.3  'C. 
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Figure  23.  Plot  of  waveform  obtained  at  steady  state 
P  =  81.2  kPa.  AT  =  394.2  °C. 
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Figure  24.  Plot  of  DFT  of  data  obtained  at  steady  state 
P  =  81.2  kPa.  AT  =  394.2  'C. 
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Figure  25.  Plot  of  waveform  obtained  at  steady  state 
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Figure  26.  Plot  of  waveform  obtained  at  steady  state 
P  =  8  0.3  kPa.  AT  =  411.4  °C. 


44 


1.2 

1 

A 

0.8 

/  \ 

E 
3 

0.6 

A 

1   \ 

5- 

0.4 

\    \ 

1             1     " 

1 

"Si 
E 
< 

0.2 

(i 
-0.2 
-0.4 

-    \ 

w 

-U.D 

0 

0.005               0.01               0.015               0.02 

0.025              0.03 

Time  (sec) 
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Figure  28.  Plot  of  DFT  of  data  obtained  at  steady  state 
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Figure  29.  Plot  of  waveform  obtained  at  steady  state 
P  =  206.5  kPa.  AT  =  395.1  °C. 
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Figure  30.  Plot  of  DFT  of  data  obtained  at  steady  state 


P  =  206.5  kPa.  AT  =  395.1  °C 

m 


46 


2 
1.5 

- 

a 

E 

D 

1 
0.5 

/- 

15 

J 
o. 
E 
< 

0 

-0.5 

-1 

-1.5 
( 

/ 

)       0.005 

0.01 

0.015 

0.02 

0.025 

0.03 

Time  (sec) 

Figure  31.  Plot  of  waveform  obtained  at  steady  state 
P  =  206.7  kPa.  AT  =  398.7  *C. 
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Figure  32.  Plot  of  DFT  of  data  obtained  at  steady  state 
P  =  206.7  kPa.  AT  =  398.7  'C. 
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Figure  33.  Plot  of  waveform  obtained  at  steady  state 
P  =  2  04.8  kPa.  AT  =  411.5  'C. 
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Figure  34.  Plot  of  DFT  of  data  obtained  at  steady  state 
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Figure  35.  Plot  of  waveform  obtained  at  steady  state 
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Figure  36.  Plot  of  DFT  of  data  obtained  at  steady  state 
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